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Abstract 

Factorization is the central ingredient in any theoretical prediction for collider experiments. 
We introduce a factorization formalism that can be applied to any desired observable, like event 
shapes or jet observables, for any number of jets and a wide range of jet algorithms in leptonic 
or hadronic collisions. This is achieved by using soft-collinear effective theory to prove the formal 
factorization of a generic fully-differential cross section in terms of a hard coefficient, and generic 
jet and soft functions. In this formalism, whether a given observable factorizes in the usual sense, 
depends on whether it is inclusive enough, so the jet functions can be calculated perturbatively. The 
factorization formula for any such observable immediately follows from our general result, including 
the precise definition of the jet and soft functions appropriate for the observable in question. As 
examples of our formalism, we work out several results in two-jet production for both e^e~ and 
pp collisions. For the latter, we also comment on how our formalism allows one to treat underlying 
events and beam remnants. 
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1. INTRODUCTION 



Factorization is the main ingredient in any theoretical prediction for colhder experiments. 
Most factorization theorems are easy to understand intuitively. For example, the most basic 
factorization theorem for the production of lepton pairs in proton-proton collisions has the 
form 

(7 = ^ (g) fi/p (g) fj/p . (1) 

id 

Here, the partonic cross section aij describes the production of the two leptons from the two 
initial state partons i and j, while the parton distribution functions fi/p and fj/p describe 
the probability of finding the partons i and j in the proton. The parton distribution functions 
and the partonic cross section both depend on the momentum fractions of the partons with 
respect to the hadrons, and the denotes the convolution in these variables. For more 
complicated processes, such as jet production, factorization formulas still exist, but are 
more complicated than for the Drell-Yan process. 

While there is usually a simple intuitive picture leading to factorization theorems like 
Eq. ([T]), many open questions cannot be answered without a much more detailed understand- 
ing of the factorization theorem. First, precise field theoretical definitions of the different 
elements in terms of matrix elements of operators are required to calculate them systemati- 
cally. Second, each of these elements depends on a renormalization scale fi, and the precise 
fi dependence cannot be obtained from the naive arguments given above. The fact that the 
final hadronic cross section is independent of fi allows one to derive renormalization group 
equations, which can be used to sum large logarithmic terms present in the perturbative 
results. Finally, it is important to understand under which circumstances the factoriza- 
tion theorems hold and when they break down. Thus, a more detailed understanding of 
factorization theorems is mandatory for a theoretical understanding of collider signatures. 

Understanding factorization has a long history, and started with the seminal work of 
Collins, Soper and Sterman (see Refs. [H,]!!, 0, 0] and references therein). It is instructive 
to remind the reader about the philosophy of these traditional factorization proofs, and to 
compare it to factorization proofs using effective fields theories, as discussed in this paper. 
While the Lagrangian of the strong interaction is given in terms of partonic degrees of 
freedom, any perturbative calculation of partonic scattering amplitudes gives rise to infrared 
divergent results. These infrared divergences are a manifestation of the well-known fact that 
at long distances the strongly interacting degrees of freedom are hadrons, not partons, and 
that the binding of partons into hadrons is a nonperturbative effect. The infrared divergences 
in the partonic results are regulated, however, if the dimension of spacetime is chosen to be 
D = A — 2e, and manifest themselves as 1/e poles with e < 0. Thus, for D 4 one can 
calculate the scattering cross section of two partons m and n, and by the same intuitive 
picture as before, one expects that the partonic scattering cross section factorizes as 

Cmn = ^ ^ ^ij ® fi/m ® fj/n ■ (2) 
hj 

In this case, fi/m [fj/n] denotes the probabihty to find the parton i [j] in the parton m [n]. 
This probability has a well-defined expression order by order in perturbation theory, and is 
infrared divergent in the limit Z) — > 4. Under the assumption that any infrared-safe (finite 
as D — 4) result is the same in the hadronic {D = 4) and partonic {D ^ 4) theories, the 
factorization of long distance and short distance physics in the hadronic theory can be proven 
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by showing to all orders in perturbation theory that aij in the partonic theory is indeed finite 
in the limit D — > 4. Traditional factorization proofs therefore use diagrammatic techniques 
to show that all infrared- divergent terms in amn are contained in the partonic distributions 
fi/m and fj/n- 

Proofs of factorization theorems in soft-collinear effective theory (SCET) [5|, Ig, l7|, l8|, on 
the other hand, use a different approach. By construction, the correct effective field theory 
reproduces the long distance dynamics of the underlying theory in a particular kinematic 
limit. SCET is constructed to reproduce the long distance physics for processes involving 
highly energetic particles, and the assumption is now that SCET is indeed the correct 
effective field theory of QCD in this particular kinematic limit. What is important is that 
this assumption can be tested in perturbation theory. 

The physics at short distances, in general, is not properly described by the dynamics of the 
effective field theory itself, but can be determined by demanding that the effective theory 
reproduces the underlying theory order by order in perturbation theory. This matching 
calculation can be performed using partonic degrees of freedom, because the effective theory 
reproduces the long distance physics of the full theory ensuring that all infrared divergences 
cancel in the matching calculation. Thus, the effective theory gives a result of the form 



a 



LD 



(3) 



where SDjj describes the short distance physics governing the scattering of two partons i 
and j into anything, while LDjj describes the long distance probability of two protons to 
give two partons i and j. The final step in the proof of the factorization formula is to show 
that 

LDij = fi/p ® fj/p (4) 

where fi/p are now matrix elements of operators defined in SCET. This is accomplished in 
SCET by exploiting the dynamics of the effective theory, as will be discussed in much more 
detail later. 

The first factorization proof in SCET was for B Dtt decays [9|, but it was soon 
reahzed that SCET can be used to reproduce factorization for simple QCD processes with 
back-to-back jets 10, 111, l^, 13|. Recently, there has been progress in studying factorization 
and the resummation of perturbative corrections for some weighted cross sections in e'^e~ — >• 
hadrons 0, E, E, 113, [li , Drell-Yan [l9|, [13, iH , deep inelastic scattering [H, H H H [25| , 
and Higgs production 2C, 2^. 

While such fully inclusive observables have proven to be very useful in capturing generic 
features of hadronic events, they are not well suited to identify specific short distance pro- 
cesses. For this reason, jet observables are usually considered, in which hadrons that are 
"close together" are collected into jets of particles. The idea is that QCD radiation will 
turn a single parton produced in a short distance process into a jet of hadrons, such that 
the total momentum of the jet can be used as a measure of the momentum of the original 
parton. Thus, jet observables can be used to directly test the underlying short distance 
process that produced the event, as long as the jets are well separated from the beam axis 
and their dependence on the underlying event is very small. 

Of course, the definition of a jet depends on the precise meaning of "close together", and 
there are many algorithms available which group the final state particles into jets [27, 2^, 2^, 



30l . l3ll . |32| . l33l . l34j . l35l . l36l . l37l . l38|. To calculate jet cross sections perturbatively, we need a 



jet algorithm that is infrared safe, such that a partonic calculation does not lead to infrared 
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divergences in D = 4 dimensions. Which algorithm is chosen by experiment ahsts is usually 
determined by practical considerations, such as speed and algorithmic robustness. 

In this paper, we develop a factorization formalism that can be applied to any desired 
observable, like event shapes or jet observables. In particular, for iV-jet production in 
hadronic collisions, we show that the cross section factorizes into a hard function, d'ij^ki...kMi 
describing the underlying partonic process to produce N partons, convoluted with N jet 
functions, J^^, a soft function, S, and parton distribution functions, fi,j/p, 

cr = ^ crij^ki...kN 'S) Jki® ■ ■ ■ 'S) JkN ® S fi/p fj/p . (5) 

The jet functions Jk^ are the final-state analog of the parton distribution functions. They 
describe how the final partons from the hard interaction evolve into the observed jets, and 
contain all dependence on the actual jet algorithm. The soft function is a nonperturbative 
object, which describes, for example, how color is rearranged to allow the colored partons to 
form color-singlet jets. The effective theory allows to give precise field-theoretic definitions 
of all objects in this factorization formula. Whether the differential cross section in some 
observable factorizes in the usual sense depends on whether the observable is dominated by 
factorizable kinematic configurations and whether it is inclusive enough to allow perturbative 
calculations of the jet functions. 

In Sec. 121 we first define a generic differential cross section written in terms of functional 
derivatives, and then show in Sec. [3] how observables are constructed from this generic cross 
section. In Sec. HI we derive a factorization formula for this differential cross section, which 
relates it to convolutions over generic building blocks. To focus on the overall structure of 
the result, this derivation will be schematic in the sense that we will ignore the explicit color 
and spin structure of the underlying interaction. As a first example, in Sec. 0, we apply 
our results to the production of two-jet events in e~^e~ collisions, including all color and 
spin information. We reproduce the known results for event shape and hemisphere mass 
distributions, and obtain results for observables based on cone jet algorithms. In Sec. O we 
apply our formalism to hadronic collisions. We derive the factorization formula for two-jet 
production in pp collisions, focussing on the simplest subprocess qq' qq' . Our conclusions 
and outlook are presented in Sec. [71 

2. ENERGY AND THREE-MOMENTUM CONFIGURATION OF AN EVENT 

The differential cross section in any observable O is given by 

^ = i E I (^1 21^) I' ^'^^y^'iPi - Px) S[0 - foiX)] , (6) 
P^ X 

where |/) denotes the initial state containing two particles with total momentum pi = 
{Ecm,0), \X) denotes an arbitrary final state with total momentum px, and the sum over 
X includes a sum over states, as well as all final-state phase-space integrations. Finally, Q 
stands for the relevant operator responsible for the underlying hard interaction. 

The function fo{X) computes the value of the observable for a given final state X, 
and in general depends on the four-momenta of all particles in X. The four- momentum 
configuration of X can be described by its energy configuration ujx{^) and three- momentum 
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configuration kx(f^)- If X lias n particles with four-momenta pi = {Ei, pj), we have 



1=1 1=1 

where fli = f^(pi) is the direction of the three-momentum p^. More generally, we can think 
of uj{Q) and k(f2) as the distribution of energy and three-momentum over the solid angle Q, 
as measured experimentally. 

To integrate over u and k, we define a functional integration measure as usual by dis- 
cretization. We divide Q into bins {fifc}, and define the set of discrete variables {ujk} and 
{kfc} as the integrals of uj{Q) and k(fi) over the bins {i^fc}/ 

Uk= f dQuj{n), kk= f dr]k(fi). (8) 

Then 

Vuj{n) = Vuj{n) e[uj{n)] = JJ dc^^ e{ujk) , 

k 

j3i, 11, |2 ^ I], I 

vk{Q) ^ vk{Q) s[Q{k{Q)) - = n ^ s{n{k,) -n,) = l[ . (9) 

The ^-functional in T>uj{Q) restricts uj{Q) to be positive, while the 5-functional in T>k(fl) 
restricts k{Q) to point into the direction Q. 

The integration measure in Eq. iQ still includes unphysical configurations. To only allow 
physical configurations, we have to include a mass-shell condition. Taking a fixed invariant- 
mass distribution as boundary condition, we get 



Vuj{n) vk{n) = / Vu{n) vk{n)5[uj{ny - k{ny - fi{nf] 

fi{n) J 

= [U^J-^^^^5icol-K\'-f^l)9M, (10) 

k ^ ' 

where fik = Jn^^^f^i^)- This fixes the direction of all particles, but could include different 
final states X, as long as they have the same invariant- mass distribution /x(r2). If we instead 
restrict the integration to a state X, having n particles with masses mj, we recover the 



This is slightly different from the usual definition of V(l){x) for some field (j){x), where the discrete variables 
4'k — 4'{xk) are taken as the values of (j) at the points xu- The difference is an irrelevant overall constant. In 
our case, taking the integrals is the more natural choice and makes the connection to the usual phase-space 
integration simpler. 
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standard n-body phase space for X, 



ix 



^^-^ J x(nu...,n„) 



i=l 



k=l 



|kfcpd|kfc 
(27r)3 



■6{ul - |kfc|^ - ml) e{ujk) 



n 



5{vl - ml) 9ip^) ^ d^x. 



(11) 



On the right-hand side of the first hne, we first integrate uj{VL) and k(i7) with the boundary 
condition that there are exactly n particles with masses mj moving in the directions 
denoted as X(f2i, . . . , f2„), which is then integrated over the particle's positions Qj. In the 
second line, in the discretization only the integrals over those n bins survive that happen to 
contain a particle. Together with the f2j integrations, this reduces to the standard n-body 
phase space for X. In the following, we will mostly drop the dependence of uj(VL) and k(f2) 
on n, but one should always keep in mind that uj and k are functions of VL. We will still use 
square brackets to denote functionals f[uj] and /[k]. 

Returning to Eq. ([6]), we now assume that fo{^) does not depend on any internal quan- 
tum numbers of X, but only on the four-momenta of all particles in X.^ In this case, fo 
can be written as a functional of the energy and three-momentum configurations. 



fo = fo[^^, k] with fo{X) = fo[uJx, k 



X 



(12) 



We now define an energy- momentum flow operator = S^{fl), whose eigenvalues are the 
energy and three-momentum configurations of the state \X) in Eq. ([7]), 



(13) 



The energy flow operator £^(£i) has been used previously, for example to study two-jet event 



shape distributions [IJ, ll7|, |39|, |40|, |4l|, |42| and jet energy-flow correlations [43|, |4J, |45|, |46 

In terms of the energy-momentum tensor 



dC 



(14) 



we can write E^iVi) as [13, 39l 

e^'iSi) 



i?— >oo 



lim i?2 / dtn^T'^'(t,i?n). 



(15) 



Here, n = n(f2) is the unit three-vector pointing in the direction identifled by fi. Therefore, 
£^^(n) measures the total four-momentum arriving over time at inflnity in the direction 
An expression for £^°(n) similar to Eq. (fT5i) in terms of an integral over R for t ^ oo 



^ Note that most information about internal quantum numbers, such as the number of 6-jets, is obtained 
from four-momentum information alone. 
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was derived in Refs. 4^, |47|. An explicit proof of Eq. (fTSll for S^{Q) for scalars and Dirac 



fermions can be found in Ref. 17 



Using Eqs. (fT2l) and (fT3!l . we can write Eq. ([6]) as 

= ^JvujVkj d^x (/|Qt(x) <5[u; - (5[k - £]Q(0)|/) 5(0 - /o[u;, k]) 

VLuVk^^5iO-fo[to,k]). (16) 

In the second line, we shift to position x, turning the momentum conservation into 
an integral over x, and rewrite ux and kx in terms of S'^. This removes any explicit 
dependence on the final state X, allowing us to perform the sum over all final states states 

Exl^)(^l = i- 

For the rest of this paper, we will assume for simplicity that X only contains massless 
particles. The extension to the general case is straightforward. In this case, fo = fo[^] only 
depends on uj, and we can integrate over k to find 

^ = lvj£siO-fol.]), (17) 



where the generic differential cross section 6cr/6uj is defined as 

6a _ 1 
5w ~ 2^ 



d^x {I\Q\x) 6[u; - £'] Q(0)|/) 



= wJwr " ^^^^'^^ ■ ^^^^ 

In the second line, we have written the result in terms of Q{p)-, the Fourier transform of Q{x). 
(To simplify the notation we use the same symbol for operators in position and momentum 
space and simply distinguish them by their arguments.) 

Equation f[T7l) can be regarded as the master formula of our formalism, and its ingredients 
are the subject of the following sections. We first discuss the functional fo[^] ^^e next 
section and then the factorization of 5a/5uj in Sec. HI Then in Sees. [5] and [6] we show how 
to combine these two elements to obtain a factorized form of dcr/dO for specific processes 
and observables. 



3. CONSTRUCTING OBSERVABLES 

The form of the functional fo[^] depends on the observable under consideration, and in 
this section we will give a few examples of how to construct fo [oj] for specific observables. 
To get used to our notation, we start with the simple example of the total four-momentum 
of the final state. Next, we consider event shapes, which are fully inclusive observables. 
Finally, we discuss jet observables, which are less inclusive and defined with respect to a 
specific jet algorithm. 
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A. Total Four-Momentum 



The total energy and three-momentum of the state X are 

n „ n „ 

Ex = Y,E,= dnooxi^), px = J2pi= dnn{n)ujx{n), (19) 
1=1 ^ 1=1 

where we used that for massless particles k{fl) = n{Q) uj{Q). Hence, we define 

pA*[cj] = j dnn^'in) uj{n) , (20) 

where nf'lQ) = (l,n(fi)). From Eqs. ([IT]) and ([H]) we get 

Jvuj^£ 6\P -P[co]) = ^Jd'x {I\Q\x) 6\P - P) Q(0) |/) . (21) 



d^P 

In the second step we performed the integration over oj and used that (for massless fields) 
P^[£^] yields the momentum operator P^ = (P°, P),^ 

P^'[£0] = j dnS^in) = hm j dxT'^°(t,x) = Pf". (23) 
B. Event Shapes 

Event shapes are defined with respect to the thrust axis of an event. Given a final state 
X, one first calculates the thrust axis t = t(X), which is defined as the unit three- vector t 
that maximizes the sum 

n 

E|t-P^l' (24) 

i=l 

which runs over all particles in X. Given t, one then calculates the observable of interest. 
A generic class of event shapes can be written as 

1 " 

MX) = J2g^(r^,(p,)) |p^(p,)| , (25) 



i=l 



where the rapidity m and transverse momentum are measured with resp ect to t. For 
example, for thrust ^M, il], jet broadening [sO], and the C-parameter 51 1, the function 



Qeiv) has the form 



^ To see this explicitly, consider the current j'^{x) = T^^{x) (for fixed 

dr2£:^(f7) = lim / dt d5n-j(t,x)=/ di / dx V • x) = lim dyif{t,x.), (22) 

«^ooJo JdS(R) Jo J t^ooj 



'dS{R) 

where in the last step we used current conservation d^j'^ix) — 0. 
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The thrust axis can be obtained from the energy configuration of the final state, and can 
therefore be written as a functional t [u] . It is defined (for massless particles) as maximizing 
the integral 

jdn\t-n{n)\iu{n) , (27) 

which for to = ujx reduces to Eq. fl21|) . Using |p^| = E/ coshr^t, we can write Eq. (12^ in 
terms of oj 

fM = Jdt6{t-t[Lo])Uuj;t], (28) 

with 

f,[uj; t] = ^ [ dngM , (29) 



Ecm J cosh r]t 

where the solid angle Q is decomposed with respect to the thrust axis as = {rit,(j)t), and 
■r]t = tanh~^(cos 6*1). 



C. Jet Observables 

1. General Features of Jet Algorithms 

A jet algorithm J' acting on a final state X returns the set of momenta of all particles 
in the event, grouped together into the different subsets belonging to each jet plus a set of 
particles not belonging to any jet, which we take to be soft: 

:r(x) = {K}i,...,K}A^,KK}. (30) 

In terms of the energy configuration uj, we can write the action of the jet algorithm as 

= {<,..., (31) 

where uf^ is the part of lu corresponding to jet i, and lu^°^^ is the remaining soft part of u; 
not assigned to any jet, such that 

u; = u;f' + ... + Jp + c^-ft . (32) 

We formally split the action of the jet algorithm into two distinct steps. We first define a 
quantity j that contains all global information about u required to construct the individual 
jet configurations Ui from u. For example, j contains the total number of jets and the 
direction of each jet. That is, j is analogous to the thrust axis in the case of event shapes. 
Hence, a jet algorithm J' provides a functional which returns the required information 

j for a given u. Second, we define functionals J7i[c<j; j] that project out the part of u belonging 
to the i-th jet. That is, for j = by definition 

Jduj;jj[u;]] = cuf , Js[uj;ij[uj]] = . (33) 

We stress that J7i[ci;; j] only encodes the actual projection, which is completely specified by 
the specifics of the jet algorithm and the information provided by j. Thus, for a given j, J7i 
can be applied to any uj. For example, by definition, J7i satisfies the consistency conditions 



9 



For simplicity, we will keep the dependence on J implicit from now on and simply write 



2. Construction of Jet Observables 

We can now write a generic jet observable as 



fo[uj] = j dj5(j - /o[^; j] 



where, in general, /©[cj; j] has the form 



(35) 



(36) 



and M is the total number of jets given by j. Most jet observables only depend on the total 
four-momentum of each jet. In this case, using Eq. fl20|) ^. we have 



i=i ^ 



9o{Pi 



P. 



N) 



(37) 



The function go{Pi, • • • , Pj\f) computes the observable of interest from the given jet momenta 
Pj. It is analogous to the function Qeiv) for event shapes. Some simple examples would be 
the total number of jets or the invariant mass of two jets. 



1, • • • , Pm 



■Af. 



gmJPl,■■■,PM) = {P^ + PJ) 



(38) 



Similar to Eq. fl37j) . one can easily define observables depending on additional information 
about the individual jets, for example, one can imagine observables which depend on a 
weighted integral of the energies of all particles in a jet. 



3. Examples of Jet Algorithms 

Of course, in practice, how the action of the jet algorithm is separated into j [a;] and J'i[uj] j] 
depends on the actual algorithm, and we will briefly discuss a few examples here. Since most 
jet algorithms do not have a simple analytic expression for generic final states, it will not 
be possible to obtain j [u] analytically, either. However, this is not a limitation, because we 
can always define j[ci;] by acting with the full jet algorithm on u and only returning the 
necessary global information. The more relevant, and perhaps nontrivial, task is to figure 
out the information required in j, and to define the projections J7^; [u ; J] accordingly. 



The simplest example is probably the hemisphere jet algorithm 28, 29, 30], for which the 
number of jets is always 2, and the only relevant global information is the axis perpendicular 
to the plane separating the two hemispheres, which is usually taken to be the thrust axis t. 
Hence, we can write j as 

j = {2;t}, i[u] = {2;t[u;]}, (39) 



^ For jet algorithms, Eq. ((20)) corresponds to the so-caUed E-scherae, which defines the total jet momentum 
as the sum of the particle momenta. 
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where by convention we included the number of jets in j. The corresponding projections J\ 
and J2 return the two hemispheres defined by the thrust axis, 

J-i[cu; j] = 0(0 < < vr/2) , jTsicu; j] = uji^l) e{'K/2 < 6, < it) , (40) 

where Q = {6^, 4>t) is given with respect to t. Note that here we have i7s[ll', j] = 0. Typical 
observables constructed from these individual jets are their invariant masses, Qm^^^Pi^ P2) = 
P^2- Another class of observables is given by the event shapes in Eq. (!29l) . for which 

A less trivial example is a cone jet algorithm [27|, l35|, ISg, [371, l38| . In this case, the necessary 



global information returned by j[ci;] is the number of jets A^, and the direction jj of each jet, 

j = {Ar;ji,...,j^}. (41) 

For instance, in e'^e~ collisions one can define an A/'-jet final state as one admitting a 
minimum number Af directions jj such that the total energy outside an opening half angle 
R about each direction is less than some fraction e. For Af = 2, and with the additional 



constraint ji = — j2, this is equivalent to the original Sterman- Weinberg jet definition [27 
For a given set j, the projections are then simply 

J,[co;j]=ujin)9iR-eiJ, (42) 

where R is the cone radius and Q = {6j-, 0j-) is given with respect to jj for each i. 

As with all observables, we require jet algorithms that are infrared safe, which is not the 
case for many cone jet algorithms. An example of an infrared-safe cone jet algorithm is the 
seedless infrared-safe (SIS) cone jet algorithm (sTj . For illustration of our method, we will 



use the snowmass cone algorithm [35|, l38| as an example, however our approach can easily 
be adapted to the SIS algorithm or any other infrared-safe jet algorithm. 

For hadronic collisions, the variables (?7, 0), defined with respect to the beam axis, have 
simple transformations under boosts along the beam direction, and so are favored over {6, </>). 
In the Snowmass cone algorithm, jets are defined by cones of constant radius R in {rj, (p) 
space. When applied to massless particles, the directions jj are given by the solutions of ^ 

= JdQ ^ [n - fi(jO] 9{R - v/[^-^a)F + [0-0(j.)]^) , (43) 

where fl = {ri,(f)) is now measured with respect to the beam axis, and fi(ji) = (^(ji), 0(ji)) 
are the coordinates of the i-th jet direction. Equation fH3l) is the analog of Eq. (127|) for the 
thrust axis. The corresponding projections are 

J,[u;}] =u;{n)e{R- - ^(j.)]' + '/'(j.)]') • (44) 



Finally, kx jet algorithms [Ml, l32, l33|, l3^ also fit our general definition of algorithms. 
Although their definition contains a cut on some distance measure, the precise size and 
shape of a jet also depends on the details of the energy configuration u. Hence, there is no 
simple expression for the projections jTjct;; j] for arbitrary u and j. In principle, they are 
well defined (albeit complicated) for a fixed number of particles, in which case integrals over 
the energy configuration reduce to normal phase space integrals [see Eq. f llll) ]. In practical 
applications, it is easiest to apply the algorithm numerically. 



^ For some configurations uj, this equation can admit multiple sets of solutions j, containing a different 
number of jets. This happens when there are overlapping cones, and one has to decide whether to split 
or merge these. 
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4. FACTORIZATION OF da/5u 



In this section, we prove a factorization theorem for the generic differential distribution 
5(j/5uj defined in Eq. f|T8|) . This will allow us to separate the various scales in the problem 
and write our result in terms of convolutions over simpler functions, each of which captures 
only the physics at a certain energy scale. The factorization proof uses arguments similar 



to those used to prove factorization for event shape distributions in Ref. [17[. The central 
ingredient in addition to the usual factorization of soft and collinear degrees of freedom in 
SCET will be the use of the energy flow operator £^{VL) defined in Eq. (ITSll . 

When deriving the factorization formula, we will ignore all color and spin structure of the 
SCET operators, and denote all collinear fields by 0, regardless of whether they correspond 
to quarks, anti-quarks or gluons. This schematic notation will allow us to focus on the issues 
directly related to the proof of factorization. Of course, to obtain the full result for the cross 
section, the color and spin information has to be included, and we illustrate how this is 
achieved when discussing explicit examples in Sees. E] and El 



A. Matching QCD onto SCET in Momentum Space 

Usually, the matching of QCD onto SCET is performed in position space by expanding 
the relevant QCD operator Q{x) in terms of SCET operators 0(a;), 

Q{x) = J2 ^'^'-'^ (^K,pj(x) . (45) 

Here, C{„.}(9{„-_p.} stands for a sum over several SCET operators with the same number 
of collinear directions, each with its own Wilson coefficient. The Wilson coefficients are 
determined by taking matrix elements of both sides, and expanding the full-theory matrix 
elements in terms of the matrix elements {0{x)) evaluated in SCET. Note also that at 

this point the operators include all incoming and outgoing fields, whether they are strongly 
interacting or not. 

The operator (9{„. p-}(a;) is written in terms of (gauge-invariant) collinear SCET fields 
(pn.pix). Each field depends on a large label momentum p^^ = p~n^/2 + with = 
and ~ (^{^P~)i and has a residual x dependence corresponding to a residual momentum 
k ~ (9(A^p~), so the total momentum of the field is p = p + k. Thus, one can think of the 
fields (f)n,p{x) as being written in label momentum space and residual position space. With 
this interpretation, the sum over all labels rti and Pi in Eq. fHS]) corresponds to taking the 
remaining label Fourier transform to convert the right-hand side to full position space. 

This separation into discrete label and continuous residual components is conceptually 
convenient when formulating the effective theory, and means that phase space is divided up 
as 

U^P = Y^ I d^k . (46) 

The concrete choice of the discrete labels n and p is determined by the external momenta. 
As is well-known, this choice is arbitrary at subleading order in A, which can be exploited 



to derive constraints from reparametrization invariance [5^, 153 . 

However, in practical applications, especially with more than one collinear direction, the 
label choice can easily get obscured during the calculation. One example is four-momentum 
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conservation for two back-to-back jets with collinear momenta pi = pi + ki and p2 = P2 + ^2, 



where 5^ ~ 0{ki^2)- Both equations are correct and correspond to different choices of the 
label momenta. Using the first equality, as is often done, may seem somewhat ad hoc, but 
one simply makes an implicit choice of, say, p2 relative to pi. Of course, this is only justified 
if P2 was not already chosen somewhere else. Furthermore, at the end of the day, one often 
has to recombine some leftover label sums and residual integrations, e.g.. 



n 



dk^ an, (48) 



corresponding to unconstrained phase space integrations of external particles. 

With several collinear directions, keeping track of all label choices and dealing with 
leftover label sums and residual integrations quickly becomes very cumbersome. To avoid 
all of these issues, it is convenient to perform the matching entirely in momentum space, so 
Eq. (145|) becomes 

Q{p) = E (n / tSi) (2^)'^'(^^ - J2(P^ + ^«)) Cin.m}) Ok,p.}({A:.}) . (49) 

Here, Q{p) is the Fourier transform of Q{x), and 0{ni,p,}{{ki}) is written in terms of 
momentum-space SCET fields, (f)n,p{k), which are obtained by taking the remaining residual 
Fourier transform of 0„^p(x), 



0„,p(A;) = y d4xe^'^-^0„,p(x). (50) 

We can imagine that the matching is performed at fixed total momentum Pi = pi + ki of 
each field in 0^ni,pi}- We then choose the field labels directly during the matching such that 
n = p/|p| for each field. With this choice, p~ = p~ = p^ + |p|, p^ = k^ = p^ — |p|, and 
p\_ = 0. This allows us to recombine the label sums and residual integrations in Eq. (149|) 
into df^Pi integrals^ 

Ef d"^^ _ f f dp~dp~^ dQ {p~ — p'^)'^ , . 

^,7(2^ = 7(2^ = 7^2^0^ 8 ' ^ ^ 

where Q = f2(p) is the solid angle corresponding to the direction of p. We also keep the 
dependence on the labels implicit and simply write the fields in the operator in terms of 
their total momentum p, 

0(p) ^ 0„,^(A;) , Oi{p.}) ^ OK,pa({fc.}) . (52) 

Hence, the final form of the matching becomes 

Qip) = (n / ^) (2^)'^'(p - E^O ^^^p'^^ ^^^p'^^ ' (^^^ 

where here and in the following it is understood that the Wilson coefficient C only depends 
on the directions rij = Pi/|pj| and large components p~ = + Pi of the momenta Pi. 
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We suppress that, strictly speaking, the integral over p should be restricted to only include collinear 



momenta, which is equivalent to excluding the zero-bin region 5J| from the integral 
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B. Factorization Proof 



Starting from the definition of bo j bio in Eq. f|T8|) . in the first step we match QCD onto 
SCET. According to Eq. fl53|) . the matching condition takes the form 



X Oi{pa,Pb) Of{Pi,...,Pn) {'^'^yS'^(p-Pa-Pb + ^Pi] ■ 

^ 1=1 ^ 



(54) 



Here, Pa and pb are initial state coUinear momenta, and the operator Oj is responsible for 
annihilating the initial state. Similarly, pi, . . . ,Pn are N final state coUinear momenta, and 
the operator Op, defined as 



N 



Op{p,,...,PN)=l[<p\p^), (55) 

4 = 1 

is responsible for creating the final state. Equation (15^ is valid in any region of multi-body 
phase space that is dominated by jets of coUinear particles, corresponding to coUinear 
directions, that are well separated from each other and the beam axis, i.e., the initial coUinear 
directions pa,b- The dominant power corrections to Eq. (JSlj) scale like p'^/pi ■ Pj. 

The different coUinear fields in Oj and Op interact with each other only through the 
exchange of soft gluons. These interactions are eliminated to all orders in and leading 
order in the power counting by the usual field redefinition in SCET [sj, 

0n,p(x) = Yn{x) (p'-^lix) , (56) 

where Yn{x) denotes the appropriate soft Wilson line along the direction n in the color 
representation of 4'n,p{x)- For color singlet fields, Y{x) = 1. As usual, we will drop the 
superscript on the redefined fields and operators henceforth. In terms of the redefined fields, 
the matching condition in Eq. fl54l) takes the form 



Q{p)=(^fl J -^^^ CN{Pa,Pb;Pl,---,PN) J 



(2vr)4 

X Oi{pa,Pb) Of{pi, ■.■,Pn) Os{ks) {2'nf5^ip-pa-pb + K + ^Pi\ 

^ 1=1 ^ 



N 

(57) 



where the soft operator Os{ks) contains the Fourier transform of the time-ordered product 
of all soft Wilson lines. 



Os{ks)= d^xe 



N 



y;,(x)F„,(x) n^nU 



i=l 



(58) 



Having factored the operator Q(p), we next move our attention to the 5[uj — term 
in Eq. (fTSl) . After the field redefinition, the leading-order SCET Lagrangian with A^ + 2 
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collinear directions can be written as 



N 



SCET 



(59) 



where Ci only contains collinear fields in the direction n^, and Cs is the purely soft La- 
grangian. Since the energy- momentum fiow operator, defined in Eq. ffTSjl . is linear in the 
Lagrangian of the theory, we have 



N 



(60) 



where S^^^^fl) is defined analogously to Eq. (IT^ . but using the energy-momentum tensor 
obtained from the Lagrangian £j only. Thus, S^{Q) describes the energy-momentum fiow 
in the z-th collinear sector, while S^{fl) describes the remaining soft energy-momentum fiow. 
This allows us to write 



5[ 



UJ 



UJ — ijj^ 



N 

E 

i=a,b,l 



(61) 



Combining Eqs. (1611) and (!57|) with Eq. (ITSl) . and letting = {paiVbWii ■ ■ ■ iVn} denote 
a point in (2 A^)-body phase space (with d$Ar the corresponding phase space measure), 
5(j/5uj can be written as 



5. 



a 



5uj 2p] 



X / 



N 



(27r)4 (27r)^ 



N 

n 

i=a,bA 



VuJi 



N 



i=a,b,l 



X (27r)M^($^-A;,)5 



N 

E 

i=a,b,l 



(62) 



Since there are no interactions between the different collinear sectors or the soft sector in 
Eq. fl59|) . we can factorize the forward matrix element into a product of several matrix 
elements. 

First, for each final state collinear sector we get the vacuum expectation value of two 
collinear fields, with an insertion of (5[a;j — between the fields, which restricts the collinear 
energy configuration to Ui. Since the matrix element conserves four momentum, we can write 
it as 

{mP'i) '^h - ^i] 0^fe)|O) = (27r)V(p: - p,) J{pf, UJ,) , (63) 

which defines the momentum-space jet function J{pi; uJi). Integrating both sides over pi, we 
obtain the equivalent definition in position space in terms of the standard fields (pn.pix), 



Jip,,uj,) = ld'xe''^-^0\<l,rr.,pAx)S[^.-£n<l>l,pM^) 



Ol a\ 



(64) 
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Momentum conservation implies that J{p\ uo) only has support for configurations uj that 
satisfy P{uj\ = p, where P[uj] is defined in Eq. ( 120!) . One can also define a jet function 
J{uj) = fd'^pJ{p;uj), which has support for any (physically allowed) u, but we find it 
conceptually and notationally easier to keep the total momentum of u as an explicit separate 
argument. 

At leading order in the power counting, any hadron in the initial state is bound by 
collinear interactions only, and thus does not interact with the soft sector. Hence, we can 
factor out the soft matrix element, 

(0|Oj(A;;) S[uj, - Osihm = {2txY5\K ' i2n)%h - ^N) 5(a;,) , (65) 

which defines the soft function S{uJs), and we again used momentum conservation. Note 
that Os and S depend on the + 2 collinear directions with respect to which the Wilson 
lines are defined, which is hidden in our notation. The soft function S{u!s) is defined with 
support for any physical Ug, and the total soft momentum is given by kg = P[uJs]- 
The remaining initial state matrix element can be written as 

{l\0\{p'a,p',) - S'J 6[u, - S^] Or{pa,pt)\I) 

= (27r)^5^(K - Pa) {27:Y5\p', - p,) I{pa, Pb, UJa, ^fe) , (66) 

and defines the initial-state function I{pa,Pb',uJa,^^b)- In writing Eq. fl66|) we already used 
that the matrix element will factorize for the two collinear sectors, which allows us to write 
two separate momentum-conserving 6 functions. As for the jet function, we choose to keep 
the momenta explicit in the definition of /(pa,pb; tUa,^^), so its support is restricted 
by momentum conservation to P[uJa + ujb] = Pi — Pa — Pb- For e^e~ collisions, the initial- 
state function reduces to a calculable leptonic matrix element, as discussed in Sec. 15 A[ For 
hadronic collisions, it can be reduced to parton distribution functions, but also allows one 
to treat the underlying event or beam remnants, as discussed in Sec. 16 B[ 

Combining Eqs. (163|) . (165|1 . and (!66|) with Eq. (!62l) . we can perform the integrals over all 
primed momenta, and arrive at 



N 

X \CN{Pa,Pb]Pl,---,PN)\^ IVuJsSiuJs) 



/ ^^^^■^(^'-^*)) j^,^4'^^-'^^bI{Pa,Pb;UJa,UJb) 



/ ^ \ r ^ 

X (27r)V( Pa +Pb - P[uJs] - j S UJ-UJ,- ^ 

^ 1=1 ' '- i=a,b, 



i=a,b,l 



(67) 



As anticipated, the fully differential cross section 5a/5uj is given by the product of a hard 
coefficient, ICat^, A^ jet functions, J{pi;Ui), an initial-state function, I{pa,Pb',^^a,^^b), and 
a soft function, S{uJs). Note that there are no power corrections to Eq. (!67|) other than 
from higher-order SCET operators in the matching of QCD onto SCET and higher-order 
contributions to the Lagrangian, which could in principle be included systematically. One 
should keep in mind that this factorization is purely academic at this point, because all 
ingredients depend on the precise energy configuration in each sector of the theory. The 
energy configurations are obviously very different for partonic and hadronic states, and 
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therefore the functions J, /, and 5* cannot be calculated perturbatively. One should think 
of them as fully exclusive functions. 

The importance of Eq. (!67|) lies in the fact that it establishes factorization for a generic 
A^-jet like kinematic configuration. In our formalism, the question whether the cross sec- 
tion dcr/dO for a particular observable factorizes is two-fold. First, a given value of the 
observable has to be dominated by factorizable kinematic configurations. If this is the case, 
one can immediately obtain a factorized form for da/dO from Eq. fl67l) via Eq. fll7l) . This 
means that any jet observable (meaning any observable whose definition restricts it to A^-jet 
configurations) is formally factorizable. The second, and more important, question then is 
whether one is able to determine the relevant functions, J, /, and S, for a given observable. 

For sufficiently inclusive observables, the jet functions, J{pi; uJi), will be smeared enough, 
i.e., integrated over Ui with a sufficiently smooth weight function, such that we can trust their 
perturbative calculation. Similarly, the soft function, S{uJs), and (for hadronic collisions) the 
initial-state function, I{pa,Pb',uJa,^^b), have to be smeared enough (integrated over and 
u!a,b) to reduce to nonperturbative functions that are universal between different processes. 
For such observables, one obtains a factorization formula in the more traditional sense, 
which allows for the perturbative calculation of all ingredients, with the exception of maybe 
a traditional soft function or initial state parton distribution functions. 

To study the structure of the factorization for a specific observable, and obtain explicit 
definitions of the relevant jet, soft, and initial state functions, it is usually required to also 
expand the kinematics of the process, because Eq. (167|) still mixes momentum components 
with different scaling in SCET. In this way, one obtains a result that formally is fully 
leading order in the power counting. As discussed above, with our choice of field labels the 
components pf = p^ ± \pi\ are defined with respect to the direction of the momentum pi 
itself, so pi± = and n, = Pi/|pi|. Since p~ ^ pf, to leading order the phase space in 



Equation (1671) together with Eqs. (168!) and (169|) provides the final factorized form of the 
fully differential cross section 6a/6uj for A^ jets, and is the main result of this paper. In 
the remaining part of the paper we will show how to use this result to understand the 
factorization properties of several observables. We will focus mostly on simple two-jet final 
states, for which the kinematics is simple enough to explicitly perform all phase space 
integrations analytically. All our examples, however, follow directly from our general result, 
and the extensions to more complicated final states should be obvious. 

5. e+e- 2 JETS 

In this section, we show how to apply the result in Eq. (l67|l to the simplest case of two- 
jet events in e^e~ collisions. The analysis simplifies considerably because of the absence 
of strongly interacting particles in the initial state, and due to the back-to-back nature of 



Eq. dSH) is 




(68) 



Furthermore, expanding the momentum conserving 6 function, we find 




(69) 
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the jets and the corresponding need for only a single collinear direction, e.g., the thrust 
axis. We first give explicit definitions of the operators O/, Cj?, and Os, including all 
relevant spin and color information, and then define all the ingredients in the factorized 
form of 5a/5uj. We then apply this generic formula to the special cases of event shape 
observables in the limit e — >■ 1 and to hemisphere jet masses, whose factorization is well 
understood [i^, 41, 42, 55, 5^. Factorization for the former was proven using SCET in 
Refs. 14, 17] and for the latter in Ref. [HI, and we show how to reproduce these results. We 
then consider the factorization of generic observables defined for cone jet algorithms, and 
obtain the definition of the relevant cone jet functions and cone soft function. In SCET, 
cone jets were previously discussed in Refs. 11, [IT 



57l | using Sterman- Weinberg cones. 



A. Generic Expression 



For e^e 2 jets, including the full spin and color information, the three SCET operators 
entering the matching in Eq. fl57|) are'' 



o'j.',{pi,P2) = r{PiKx'{-P2), 



(70) 



where x'^ip) = denotes a noninteracting collinear quark field of color c and charge 

eQf (where / denotes the flavor), moving in the p/|p| direction. Note that we are distin- 
guishing particle and anti-particles by the sign of the momentum argument on the field. The 
soft Wilson lines along an outgoing collinear direction are^ 



PCX) 

FJ(x) = Pexp i^f.; / dsn-As{x + sn) 
L Jo 



(71) 



where P denotes path ordering. The Wilson coefficient at tree level is given by 



C2{Pa,Pb]Pl,P2) 



ie^Qf 

'^Pa-Pb 



[1 + 0{as)] 



(72) 



Since the initial state is not strongly interacting, the initial-state function in Eq. (!66ll . 
including the average over initial spins, reduces to 



I^''{Pa,Pb;i^a,UJb) 

4 ^ ./ (27r)4 (27r)4 



spins 



,e e e 



{Pbh'ei-P'a) - S'J 6[u;, - S^] e{-pa)Ye{Pb)\e'^e-) 



{27Cy6{pa - Pe+) {27Cy6{pb - Pe-) L 



(73) 



^ We only give the result for an intermediate photon here, and include the Z boson contribution later. 
® For a discussion^of the different choices of boundary conditions for in- and outgoing Wilson lines see for 
example Ref. 
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where Pe± are the momenta of the incoming leptons and 

L^" = p',-P:^ + p'^^p:- - g'"" {Pe- -Pe^) (74) 

is the well-known leptonic tensor. Note that as we do not consider any initial state radiation 
from the incoming leptons, = in Eq. fl75]) . 

Using (pipi) = XaiPi) 0(^2) = P2) in Eq- (163!) (where a, j3 are spinor indices), 
the quark and anti-quark jet functions become (with the sum over spins implicit) 

4'^(P2;CU2) = 1^(014,(^2) 5[^2-^2°] Xj(-P2)|0> =5'^''^(^)^^J7P2;C.2), 

where the spin-singlet and color-singlet jet functions are defined as 

Apuu^i) = 7^ / M-M^\^Mp'i) - xiPi)\0) , 



JiP2;u;2) = ^J §^M^\X{-P2) S[U^2 - S',] ^2X(-P2)|0> . (76) 

Here, tr denotes the trace over spin and color indices and is the number of colors. At 
lowest order in perturbation theory, we have J{p]uj) = J{p;u) = 2tt 6 {p~^) 9 (p^) 6 [uj{Q) — 

p'6{n-n{p))]. 

From Eq. (l65i) . the soft function is defined as 

SZfi^s) = ]^ (olO^'iK) - S^] Of{Q\0) , (77) 

where we made explicit the dependence of 5" on the directions ni^2 of the Wilson lines in Os, 
and the factor l/N^ is included by convention. Contracting with the trivial color structure 
of the jet functions in Eq. (1751) . we obtain the color-singlet soft function 

Sn.nA^s) = ^'"^''"sZfM = ^tr(0|T[(y4F„J(0)] 6[uj,-8',]T[{YlY^mM ■ (78) 

Since the spin structure of the jet functions in Eq. fl75l) factorizes, we can contract all 
vector and spinor indices, 

where Ecm is the total energy and 9i^2 are the angles of pi,2 with respect to the e~^e~ beam 
axis in the center-of-mass frame. Thus, combining all pieces with Eq. (!67j) . we find 



Sa _ Sjr'a' Qj f d^i ^ r/ . \ f 



Vui J{pi] UJi) I ^^—^1)^2 J{p2]UJ2)H2{pi,P2) [l - COs9i COS6I2) 



5u; Et^ J (27r)4 -7 (27r)^ 

X j VUs Sn^n2 i^s) (27r)^(5^(pe+ + Pe- " P[^^s] - Pi - P2) ^[^ - Ui - UJ2 - ^s] , (80) 
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where the hard coefficient H2{pi,p2) = 1 + 0{as) is defined by 



\C2iPe+,Pe-;Pl,P2) 



H2iPl,P2) ■ 



(81) 



Equation (IHOj) speciahzes Eq. (!67j) to generic 2-jet configurations u in e'^e" colhsions. 

Next, we expand the kinematics. Using Eq. (169|1 . the momentum conserving 6 function 
becomes 



Hp 



S \Pe+—+Pe 



Tlf. 



Pi 



5{p^ - 



2 

KV2 



n2 



P2 



Ecm) 5(C0S^1 + COS ^2) ^(01 



(82) 



where as before in the center-of-mass frame pi = Pe+ + Pe- = {Ecm,0)- The S functions 
allow us to perform the p^, and ^2 integrations in Eq. flHU]) . and imply that pi and p2 
are back-to-back, as expected for two-jet events. In particular, rii = —112, so we can write 
Pi^2 in terms of the components (p^,p~,n) as 



Pi = {pt^ ^cm, n(l])) , 



P2 



Er 



-n{n)) 



(83) 



where fl = {6, 0) = Qi describes the orientation of the momenta relative to the beam axis. 
We also write Snin2 = 'S'n(n)- Similar to the Wilson coefficient C2, the hard coefficient 
H2{pi,p2) does not depend on the small momentum components p^2- Since p^ = Ecm, 
we define H2{Ecra) = -^^2(^1,^2)- Combining everything with Eq. flHOj) . using Eq. fl68|) . and 
writing the momenta in terms of their components, we have 



6a 

6uj 



where 



H2{Ecra 
, dcTo 



2^ 



dcos^^ 
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Vus Snip.) iy^s) 5[u: - ui - u:2 



2tt 



Vuj2 J{p2,Ecm, -n{VL)]u;2) 



dcTn 



TTtt 



(84) 



(85) 



dcos^^2Ei:^^^Hl + -s^^) 
is the Born differential cross section. The exchange of a Z boson can be included by using 



dan 



ira 



dcos^ 2EI 



+ 



Nr. 



2v,VfQf , iv^,+al){v} + a}y 



Q 



n2 



+ 



AaeajQj 



ml/E^ 



2 

cm 



- ml/E, 

8 Vea^Vfaf 
'l-ml/E^^ 



;i - ml/E^^ 



cos 9 



1 + cos^ 



(86) 



where Vej and Og,/ are the standard vector and axial couplings to the Z. 

Equation fl84l) is the penultimate formula for generic observables in e~^e~ 2 jet events. 
Each of the ingredients in Eq. (18^ is a completely exclusive object that depends on the 
energy distribution of the individual partons. The details of how to integrate over the 
energy configurations to arrive at perturbative jet functions and a universal soft function 
depend on the observable in question, but since all observable independent simplifications 
have been done, a wide range of factorization theorems can now be obtained with relative 
ease. We illustrate this with a few examples in the following subsections. 
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B. Event Shapes in the Limit e — >^ 1 

Combining Eqs. ffTTl) and fl28|) . the differential cross section in some event shape e is 

^ = jdtjvuj^dit-tiuj]) 6{e - fe[u- 1]) . (87) 

For e — i> 1, the final state is dominated by two highly collimated jets, and hence, we can use 
the result for Sa/Sco in Eq. flMl) . The integration over u is trivial and sets u = uji + uj2 + uJs- 
Since Ui^2 describe coUinear energy configurations along ±n(fi), we have t[uji + + ujg] = 
n + O(A^) This allows us to integrate over t. 



dcr fdQ fdv^ fdv~^ 

— = H2{E,^) — ^VuiJipt,E,^,n{Qy,u,) / ^ I?a;2 J(p^ , ^cm, -n(fi); c^s) 
de J 2tt J Ztc J ZiT 

j Vujs Sn{n)i^^s) S{e- fe[uJi + UJ2 + ^s] n(fi)]) . (88) 



d(To 
X 

dcoso' 



From Eq. ( |29il . we see that /e[w; n] is linear in uj, from which it follows that we can write 
fe[uJi +uj2 + ujs;n] = fe[uJi; n] + fe[uj2; n] + f^lus] n]. This implies 

5{e- fe[uJi +UJ2 + UJs] n]) = j dei de2 de^ 5{e - ei - 62 - 

X 5(ei - fe[uJi; n]) 5{e2 - fe[uJ2\ n]) 5(e^ - fe[u}s\ n]) , (89) 

which separates the uj dependencies in Eq. flHHl) . We stress that this is not a requirement 
for the factorization of dcx/de, as demonstrated by Eq. flHHj) . In fact, the full event-shape 
functional fe[uj\ = /e[i^, t[cij]] is not linear in uj and does not obey a similar separation, 
because t[ci;] is not linear in uj. The crucial ingredient for the factorization is the linearity 
of the energy-momentum tensor and the resulting separation of the energy fiow operator in 
Eqs. f lBUl) and fl^ . However, without Eq. fl5^ the jet and soft functions depend on the full 
energy distributions cjj, and are therefore neither perturbatively calculable, nor universal 
enough to be extracted from data. The important point about Eq. fISI?]) is that it allows us 
to perform the uj integrations in Eq. (l88|l . and to define inclusive event-shape jet and soft 
functions 

-^(^i) = / ^^^^^i</(P]^,^cm, n;u;i)(5(ei-/e[t^i;n]), 
J zn 

/dp^ 
—^Vuj2 Jivt, Ecm, n; UJ2) 6{e2 - fe[uJ2; n]) , 
Zn 

S{es) = J Vus Sn{uJs) 6{es - fe[i^s] n]) . (90) 

With the definitions in Eqs. fl75]) and fl75]) . these are identical to the definitions given in 
Ref. For J(e2) we used that fe[uJ, n] = fe[uJ, — n] because the sign of the thrust vector is 
irrelevant. By rotational invariance, after integrating over u!i^2, the jet functions J(ei), J(e2), 
do not depend on the value of n on the right-hand side. This would not be true if the thrust 
axis n in /efcu;!!] would be different from the momentum direction n in J{p'^ , Ecm,i^', uj). 
Similarly, after integrating over u^, the soft function S{es) is independent of n, because the 
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direction of the Wilson lines in S'n(co's) coincides with the thrust axis. Thus, using Eq. (pOj) 
and integrating over Q, we obtain the final result 

— = H2{Ecm)cro / deideade, J(ei) J(e2) S'(es)(5(e - ei - 62 - e,) . (91) 



where cxo = /d cos 6* dao/d cos 6' is the total Born cross section. Equation (iQljl agrees with 
the result of Ref. 17 . 



C. Double Differential Hemisphere Mass Distribution 

Combining Eqs. f|T71) and fl37|) . the double differential hemisphere mass distribution is 
d^cr /",. /",. dV 



dMf dM| 



I d'P, I d'P, ^j^^ 5(M,2 - P^) 6{Ml - P|) , (92) 



where the cross section differential in the total momentum of each jet for the hemisphere jet 
algorithm is 

dV 



^4p^^4p^ =/dt y"Pc.^5(t-tM)m-^[:^ikt]])5V2-m2kt]]). (93) 
Here, P[uj\ is given in Eq. (1201) and Ji[uj\ t] in Eq. (l40l) . Combining these we have 

PLnA^-^] = P^JM,-^]] = j dnn^{n)uj{n)e{() <e,< 71/2) , 

n''emi2[^; t] = ^''[J2[u;; t]] = j df]n^(l]) C^(fi) 9(71/2 < Ot < Tl) . (94) 

In general, Eq. (1^5]) will receive contributions from final states containing several distinct 
collinear momenta, corresponding to SCET operators with N > 2. However, if the final 
states are restricted to the kinematic region of small hemisphere invariant masses = 
Pf ~ 0{\^El^, corresponding to two coUimated jets, the operator with N = 2 collinear 
directions gives the dominant contribution, with the corrections suppressed by powers of A. 
Thus, we can apply the result for Sa/Su in Eq. flH^ in this region. 

The integral over t can be performed as in the previous subsection, which sets t = 
n + 0{X^). As Eq. fl94l) is linear in u, we then have 

-Phemi4^; n] = Phomii[^l +(^2 + ^3;^] = Phomii[t^i; n] + Phomii[t^2; n] + Phomi^^s! . (95) 

To understand the size of Phemij[^j] ^ = J ^^id z 7^ j, we need to think about states in 
SCET in some more detail. Since the direction n labelling the collinear fields in SCET is a 
conserved quantum number, there exists a basis for the physical states which have a fixed 
value of n as well. This implies that for a given SCET state with momentum p, one has to 
identify the value of the direction n as well. Of course, to have the same final states as in 
full QCD, one needs Jt-)^*"^""" = 1^)^*"°, i.e. one has to be careful not to double count 

the physical states. Certainly, a convenient choice is to define the SCET states such that for 
every momentum p there is only a single value n. For our problem, the simplest choice is to 
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assign the label rii = (1, n) to all states with momentum in hemisphere 1, and n2 = (1, — n) 
to all states in hemisphere 2. This choice implies 

-PhemiiKvi; n] = and Phcmii[uJi;n] = pi , (96) 

where pi = P[uJi] is the total momentum of Ui, i.e., the momentum in J{pi\uji). The power 
counting of SCET implies that (ii = Phemii[<^s] ~ '^^^Irai where Ci can be interpreted as the 
total soft momentum in each hemisphere. Thus, using Eq. (!83l) we can expand 

Mf = Pf= p^ + 2p, ■ + £2 = E^^M + ni ■ + 0{\'El^) . (97) 

Since our observables Mf only depend on pf and rii ■ ii, we can do the remaining inte- 
grations in Eq. (18^ . and define the corresponding jet and soft functions 



2vri^cm J 

J (E^mPt) = r, F / J {pt, Ecm, n; UJ2) 

ZTTEcm J 



>S'hemi (^1 5 ^; 



2 



j VUs SniuJs)S{it - ni-Phemil[uJs; n]) 5(£^ - ^2 " Phemi2 [t^^; n]) . (98) 



Again, after integrating over u!i^2,s, the jet functions, J{pl) and Ji^pl), and the hemisphere 
soft function, Si^emii^t , ^t) y do not depend on n due to rotational invariance. The SCET 



hemisphere soft function has been discussed previously in Refs. pJl, Il6|, |59|, 160|. The above 
definition provides an operator definition of S^eraii^t , ^t) SCET, and is equivalent to the 
definition in Ref. jllj. Putting everything together, we obtain for the double differential 
hemisphere mass distribution 



dV 



dMf dMi 



i^2(^cm) CTo J dit dit J{Ml - E,^iX) J{Ml - E,^et) 5hemi(^^ 4) , (99) 



which agrees with the massless limit of the result in Ref. [15 



D. Two-Jet Cone Algorithms 

As the last example in this section, we consider the cross section for two-jet final states 
obtained from an infrared-safe cone jet algorithm. Since the discussion follows closely that of 
the previous two subsection, we keep it short, mainly highlighting the differences. Combining 
Eqs. (inD and ([37]) with A/" = 2 we have 

^ = y djidja j Vuj^5{ji -ji[uj])5{i2-h[^]) 

X y"d^Pi5^(Pi-Peone[^;jl]) J d^P2 S\P2 - Pcono[uJ; h]) ^{0 - go{Pl, P2)) , (100) 

where }i[uj] denotes the i-th jet direction returned by j[a;], and the functionals for the total 
jet momenta are now defined for example using Eq. fH2l) 



^c'onekj.] = = ldQn^iQ)uj{Q)e{R-ei^). (101) 
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As before, Eq. (llOOp receives in general contributions from operators with N >2. However, 
if the final state is restricted to two jets with small invariant masses, Mf = Pf ~ 0{\^E'^^, 
the result for 5a /5uj in Eq. (18^ for N = 2 gives the dominant contribution, with corrections 
suppressed by powers of A. The restriction on the kinematics of the final state is now 
provided by the jet algorithm, or by the combination of jet algorithm and observable. 

For a good jet algorithm, the result of j[a;] should not depend on Ug up to power correc- 
tions. This is equivalent to the requirement that the jet algorithm should not be infrared 
sensitive. Furthermore, since uji describes a coUinear energy configuration along rij, by a 
similar argument as in the case of thrust, up to power corrections, the direction of the jets 
is aligned with the direction of the collinear fields. Therefore, 



-^i[uji+uj2 + ujs] = ni + C(A'=) 



(102) 



The power of k depends on the details of the algorithm, e.g. for the hemisphere jet algorithm, 
where j is the thrust axis, we had k = 2. 

To define the states in SCET, we assign the label Ui to states with momentum lying in 
the i-th cone, so there is again no overlap between states with the same momentum but 
different n inside the cones. The precise definition of states with momentum outside any of 
the cones is not important at this point. With this definition, Pcone[^jj^i', Hj] = 0, and since 
Eq. (11011) is linear in a;, we have 



(103) 



where Qi — /'cone['-'^ii Hi] &nd ii — -Pcono [uJs] Hi] are the total collinear and soft momentum in 
each cone. Equation fll02p implies that Qi+ii are aligned along n, up to power corrections. In 
addition, note that qf = qf{R) is a function of the cone size R (and the used jet algorithm). 
For R = IT the cones become hemispheres, and thus qfiir) = pf, while at lowest order 
in perturbation theory, qf{R) = pf — 0{as). Thus, for large enough R, qf /pf ~ 1 with 
the corrections calculable in perturbation theory. (Generically, we expect the perturbative 
corrections to contain logarithms of vr/i?. Similar phase space logarithms have been studied 
for the case of Sterman- Weinberg jets in Ref. |57|.) Hence, as qf obeys the same power 
counting as pf for reasonable i?, any observable that does not vanish at leading order in the 
SCET power counting can be written as 



9o{Pi.P2)=9o{qt 



'{,qT,(i2 + 



"2-1 ^2 



n 



(104) 



where, m is not necessarily the same as k and also depends on the observable. 

Since Eq. (11041) only depends on qf, the result for da/dO can be expressed in terms of 
the cone jet and soft functions 

/dp"*" 
-^VuJiJ{pf,Ecra,'n;uJi)5{qf - nrPcone[c^i; T^])5{qi - n2-Pcone[c^i; n]) , 

"dp^ 



Jcone{,q2 i 0.2 



2tx 



[uJ2]VL])5{q2 -n2 ^2; n]) 



5'cone(^^, ^t)= I T^^s Sr,{uj,)5{lf - ^ 1 " Pcone [t^s ; n]) 5^(£^ - 7^2 • Pcone [t^s ; -n]) 



(105) 



where as before ni 2 
n. Combining Eq. 



(l,±n), and the functions on the left-hand side do not depend on 
with Eq. (11001) and using the above definitions, we obtain the final 
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result for the factorized differential cross section 



= H2{Ecm) [ ^ / ^^1 Jconciqt,C[l ) / dgj dq^ Jconc(gJ,g2 



dO J 2-K dco^e 

X / d£+ d4 S,onc(^^ 4) '^(O - go{qt + gr, ^2+ + 4, 92^, n(r]))) . (106) 



To our knowledge, factorization forjet distributions has not received much attention in the 
literature (however, see Refs. 57|, l6l|), and this is the first time any factorization theorem 
for jet observables based on jet algorithms has been proven in the framework of SCET. 

For many observables, such as the transverse momentum distribution, the dependence 
on the soft momenta and the small components ql is power suppressed, which allows us to 
integrate over these to obtain 

^ = /f2(^cm) / ^ / ^^l" '^cono(gr) / dq^ Jcone(g^) 5'cone 5{0 - ^©(C 9^, , 

dO J 111 dcost^ J J 

(107) 

where / d/"'"d/~ 5'cone(^''', ^~) = 5'cone is perturbatively calculable up to small power correc- 
tions and we defined 

^cone(gr) = j Ml Jconciqt, Ql) , Jconc{q2) = j d^^ J cone {qt , '?2') • (108) 



6. TOWARDS pp^2 JETS 



In the previous section we have focused on two-jet production in e'^e~ collisions. In 
this section, we extend these results to include hadrons in the initial state. Jet production 
in hadronic collisions is in several ways more complicated than for e~^e~ collisions. First, 
there are several different partonic processes contributing to pp ^ 2 jets. Second, the 
operators describing the short distance process now contain strongly interacting particles 
for both initial and final states, giving rise to a more involved color and Dirac structure. 
Finally, there are several additional matrix elements required to describe the long distance 
physics. These are the parton distribution function describing how the initial state partons 
are distributed inside the incoming proton, as well as new soft functions. 

In this paper, we will only consider the simplest partonic process qq' — *• qq', and work 
only to tree level in the matching from QCD onto SCET. This simplifies the discussion 
dramatically, since only a single operator contributes at this order. Furthermore, due to the 
absence of gluons in the initial or final state, the only additional nonperturbative ingredients 
are the parton distributions of finding a quark inside the proton and the new soft function. 
The complete analysis of pp ^ 2 jets is considerably more involved and will be discussed 
elsewhere j62| . 

A. Matching onto SCET at Tree Level 

At tree level, only a single operator is required in SCET to describe the partonic process 
qq' ^ qq' 1 schematically 

C>{Pa,Pb;Pl,P2) = Ci{pa,Pb;Pl,P2) Ol{pa,Pb) Cf(Pi,P2) Osih) , (109) 



25 



where we define the Wilson coefficient C4 to contain all the kinematic and Dirac factors. 
The operators O7, Op, and Os are defined as 

oUp^,p.) = rM)xiiP2), 



rib J 



X 



;iio) 



where subscripts denote spinor and superscripts color indices. The Wilson lines for the 
outgoing fields are defined as in Eq. (1711) . while for the incoming fields they are 



Yn{x) = P exp 
The Wilson coefficient is given by 



ds n ■ As(x + sn) 



;iii) 



Ci{Pa,Ph]Pl,P2) = -J- {l'')ya {lii)sfS 



(112) 



where we stress again that we are only working to tree level in the matching. The variable 
t is one of the usual Mandelstam variables defined in terms of the partonic momenta 



s= {Pa + Pbf , i = {pa - pif 



U= {Pa- P2) 



;ii3) 



B. New Nonperturbative Matrix Elements 

There are two sources of additional matrix elements which cannot be calculated perturba- 
tively. First, the operator Oj now includes strongly interacting degrees of freedom, and the 
matrix elements involving the initial state protons are no longer calculable. Second, the soft 
operator contains four Wilson lines, rather than just two as for e^e~ collisions. This implies 
that a new soft function is required. In this section we define all required nonperturbative 
matrix elements needed for the process pp ^ 2 jets via the partonic process qq' — > qq'. 

Since the initial state hadrons are moving along different light cones, they are described 
by two sets of SCET Lagrangians which do not interact with each another. Therefore, the 
physics of the two initial states completely factorizes, in the same way as the final state 
partons in different directions factorize from one another, and we can write |J) = \Pa) \Pb) 
and Oj{pa,Pb) = 0'}{pa)0^j{pb), such that we can factorize the initial state matrix element 
as 

{l\0\{p',,pi) Slooa - S'J S[uj, - S^] Oj{pa,p,)\I) 

= {PalO'^^p'J 6[uJa - £'a] OnPa)\Pa) {Pb\0'j\pi) 5[uj, - O'MlP,) . (114) 

For the case considered here, the operators Oj and contain just a single quark field, 
{^i)a = XaiPa) and (0/)|^ = XpiPb)- The resulting matrix elements define the parton 
distribution functions to find the quarks q and q' in the initial protons Pa,b- 

I ^[^'^ - Xl{Pa)\Pa) = ^ S^'^ (^)^^^ fMPa, ^a) , (US) 
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and similarly for fqi/p{pb]^b)- (Note that since we are distinguishing particles and antipar- 
ticles by the sign of their momentum, there is no anti-quark distribution on the right-hand 
side.) Combining these results gives the initial state function 

li^%{Pa,Pb,^a.U;b) = ^,5''5'''{^)^ (116) 

In most cases of experimental interest, the observable is independent of Ua^b and the plus- 
and transverse components of the coUinear momentum, which means we can integrate over 
these to obtain the standard parton distribution function [l^, 13, 63l 



T^^a fq/p{Pa; ^a) 5(Pa " XaEcm) , (117) 

while everywhere else up to power corrections we can use 

Pa = XaEcm y , Pb = XbE^^ y , (118) 

with Ua = (1, np^) and nb = {l,np^) now aligned along the direction of the incoming protons. 

While the dependence on u in our generalized distributions fq/i{p; uj) is not of relevance 
for most processes of interest, it describes the energy configuration of the remnant of the 
proton after the hard scattering. Thus, this matrix element provides a field-theoretical 
definition of the beam remnant. In particular, this means that the effect of the beam 
remnant is properly taken into account in our factorization proof in Sec. HI In principle, 
operators Oj with more than one coUinear field in the directions and Ub can be included as 
well, and would describe multiple scatterings of partons originating from the initial protons. 
These additional hard scatterings give rise to what is usually referred to as the underlying 



event [6J, |65|. Thus, these effects are also taken into account in our factorization proof. 
Moreover, our formalism provides a field-theoretic basis to study the underlying event. The 
details are left for future work. 

Since we only work to tree level in the matching from QCD to SCET, there is only a 
single soft function required for the process qq' — > qq'. After contracting with the color 
structures of the initial state function and the q and q' quark jet functions Jy^(pi;i^i) and 

j/,/(p2; ^2) (defined in the first line of Eq. (1751) ). we obtain 

X T[{YlT^Y^J\Q) {YlT^Y^:)f\0)\ |0> . (119) 



2 



N,Cf 



27 



C. Generic Expression 

Combining Eqs. ([67]), ([75]), ffTT2|) . ffTTSD and ffTTOD . the qq' qq' contribution to 6a/ 600 
for 2-jet production can be written as 

Sa 1 / TT /" dVi ^ \ 1 



6uj 2EI 



/d^Pi \ ]^ 
.-a,o,i,z ^ J c 

X (27r)^5^(p„ + - P[t^s] - Pi - P2) '^['^ - tUa - t^fe - t^i - - t^s] , (120) 

where (at tree level in the matching) 

1 1 2^^ 

Hi{Pa,Pb;Pi,P2) = -tT[ijiajf,ijli-fu]-;tT^[ibl^i2Y] = ijia-ni, Tli • + ■ nft-ni). 

4 4 

(121) 

As discussed before, most observables are independent of the energy configurations uja 
and tUb, i.e., fo[^a + ojb + uJi + UJ2 + = fo[<^i + 1^2 + i^s]- Therefore we can drop these 
beam remnant configuration in the S functional for u and integrate over them in the parton 
distribution functions. Furthermore, inserting 

1 = E^^ / dxa / dxbS{p~ - XaE^^) 6{p~ - XbE^^) , (122) 
Jo Jo 

and using Eqs. flll7p and (11181) . the expression for 5a/5uj becomes 

^ = 1^/ dx„da;b/g/p(xj/g//p(xfe) j ^^yj P^i J(pi; o^i) j ^^Vuj2 J{p2;^2) 

X Hi[EcmXa^,EcmXb^]Pl,P2^ jvUsSnantmnii^s) 

X (27r)^5^ (EcmXa Y + EcmXb y - P[Lds] - Pi - P2^ S[uJ - Us - Ui - UJ2] , (123) 

where Ua^b = (l^^iPab) ^ire now aligned with the directions of the incoming protons. As in 
Sec. [5l this can be simplified further by expanding the kinematics. After some algebra, we 
obtain 

6a f dflr 



^ [ dxadxb [^^VuJiJ{pt,p^,ni;uJi) [ Vuj2 J (pt , P2 , '^2; 
ZTC Jq J 2-K J Ztc 



UJ2) 



daQ 

X fq/p{Xa) fq'/p{Xb) -. V / "^^^ '^n^ribmna (^s) S[uJ - Us - UJl - UJ2] , (124) 



d cos 9p 

where the angular integral is defined in the center-of-mass frame of the partonic collision. 
Both the large p~ components and the directions iij are functions of the partonic center-of- 
mass angular variables, Qp, and the energy fractions of the incoming partons, Xa,b- They are 
defined by 

Ecm 

Pi {Qp,Xa,Xb) = ^^[a;a(l + cos6'p)+a;fe(l- COS 6'p)] , 

P2 {Qp,Xa,Xb) = [Xa{l - COsOp) + Xb{l + COS 6p)] , 

ni{Qp, Xa, Xb) = n{Qi) , (125) 
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where fli = {6i, (pi) are given by 

Xail + cos 9p) — Xb(l — cos 9p 



cos 9i 



Xa(l + COsOp) + — COsdp] 



a 3^a(l - COsOp) - Xb{l + cosOp) 



— COsOp) + + COsOpj 

Finally, the differential cross section dcro/dcos6'p is given by 

ddo TTalCp 1 4 + (1 + cos6'p)^ 



(127) 



d cos 6'p 2Nc XaXbE^^ {1 ~ cosOp)"^ ' 
which agrees with the well known expression in terms of the Mandelstam variables s, t, u 

d(To Tca^Cp + 



D. Jet Observables 



(128) 



As an example how to use Eq. fll24p , we derive a factorized cross section for infrared-safe 
cone jet observables, which was also studied in Ref. j6l|. The required steps are very similar 
to the derivation given in Sec. 15 Dt and we only highlight the differences that arise from 
having protons in the initial state. First, 5a/5uo depends on the parton distribution functions 
fq/p{.Xa) and fq'/p{xb). Second, while for e~^e~ collisions one often uses the variables 9 and 
to denote the direction of jets, in pp collisions it is more appropriate to use the rapidity 
instead of the angle 9, due to the easier transformation properties under boosts along the 
beam direction. This gives cone jet functions Jcone that have exactly the same form as 
in Eq. fll24p but use the corresponding cone projections in place of Eq. (H2l) to define the 
functional Pcone[i^; ji] in Eq. fllOll) . The final difference is that the cone soft function now 
explicitly depends on the orientation of the directions iij relative to the beam axis, since 
it contains Wilson lines in both the directions of the incoming protons and the outgoing 
jets. In particular, this implies that the nonperturbative physics described by this cone soft 
function depends on the rapidities of the outgoing jets. The experimental determination 
of the soft function is thus considerably more difficult for hadronic collisions than for e^e~ 
collisions. 

Keeping in mind these differences, we can follow the same steps as in Sec. 15 DI to obtain 
the factorization formula for generic two-jet observable using cone jets: 

T^= / dXa dXb fq/p{Xa) fg'/p{Xb) , 

dO J zn Jq a cos 9p 



X j dq{ dq^ Jcondqt^Qi ) J dqj dq2 Jcone (g^,g2 ) 

X I d£t dit SZZi^t, It) 5{0 - 9o{qt + gr, qt + it, , ni, n^)) , (129) 



where the cone soft function is now defined as 



SZZ^it^^t) = JvUJsSn^n,n,n,{iOs)6{£i -ni-P,one[uJs]n,])6{et -n2-Pcone[uJs-,n2]). (130) 
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As for e~^e~, many jet observables only depend on the large momentum components 
and the direction of the jets. In this case, we can perform the integrals over and . 
Integrating over if we define 

Vcl45^°?n.(^r,4) = 5— , (131) 

which is now perturbatively calculable up to small power corrections. We obtain 

^ " /"^ / dXadxb/5/p(xa) U/p{xb) ^^^^ j dq^ J{q^) J dq^ 7(^2") 

X 5^°- 6{0 - go{q^, q2, n^, n^)) , (132) 
where the jet functions integrated over pf are defined as in Eq. (11081) . 



7. CONCLUSIONS AND OUTLOOK 

We have developed a new formalism for obtaining factorization theorems for almost any 
observable of interest at high energy colliders. We argued that any observable differential 
cross section can be written in terms of two building blocks, a fully differential cross section 
describing the energy and momentum distribution of a given event, together with the restric- 
tion of how to obtain the desired observable from this distribution. For events containing 
only massless particles in the final state, the only information required to define observables 
O is the energy configuration uj of the event, and we therefore focused on the cross section 
fully differential in u, which we denoted as 6a/6uj. By integrating this energy distribution 
with an appropriate functional /o[i^], the differential cross section da/dO in any observable 
O can be obtained. 

Our main result is the proof of factorization for the fully differential cross section, 6a/6uj, 
using soft-collinear effective theory. It relies on the fact that 6a/ 600 can be written directly in 
terms of a matrix element of well-defined operators in SCET using the energy flow operator. 
The linearity of the energy flow operator allowed us to factorize 6a/6uj into simpler building 
blocks, each of which is defined by matrix elements of operators in the effective theory and 
contains a single scale allowing for a systematic program of logarithmic resummation. After 
the factorized form of Sa/Su for a given process is determined once and for all, it can be used 
to study the factorization properties of specific observables. The question of whether a given 
differential cross section da/dO factorizes in the traditional sense depends on whether the 
form of fo['^] is such that it smears the individual matrix elements in 6a/ 600 into objects that 
can be either calculated perturbatively or determined experimentally from other processes. 

Using our formalism, we were able to directly study the factorization properties of the 
fully differential cross section, independent from the observable-specific functional /o[i^]- 
While the question of whether the differential cross section in a given observable factorizes 
in the traditional sense still needs to be asked on an observable-by-observable basis, this 
disentanglement demonstrates to what length the steps taken in factorization proofs are 
observable independent. It turns out that it is the observable independent analysis that 
requires most of the calculational work. The fact that we can study factorization on an 
observable independent level could potentially be relevant for Monte Carlo event generation. 
It should be possible to make a connection between our factorized result for 6a / 6uj for generic 
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A^-jet production and the A^-body partonic calculations that were introduced in Refs. [66, 63 
as input for an event generation framework. If so, our results could be used to provide 
improved theoretical inputs for event generation. However, more work in this direction is 
needed. 

To demonstrate the simplicity with which factorization formulas for specific observables 
can be obtained from the factorized result for Sa/Suj, we have applied our results to several 
simple observables in e~^e~ — >■ 2 jets. We first reproduced the known results for event shape 
and hemisphere mass distributions, and then obtained factorization formulas for generic 
observables defined in terms of the total jet momenta obtained from cone jet algorithms, 
which so far have not been studied in SCET. We have also explored some of the issues 
arising in jet production in hadronic collisions by studying the partonic subprocess qq' —>■ qq' 
using tree level matching from QCD onto SCET. In particular, we showed that the more 
complicated structure requires a soft function that is more complicated from the case of 
e~^e~ scattering. We also showed how parton distribution functions arise in our formalism, 
and commented on how it could be used to study beam remnants and underlying events. 

It should be clear from these examples how our generic N-jet formalism can be applied to 
the study of observables in more complicated processes, such as processes with heavy vector 
bosons and more than two jets in the final state, which are crucial for many measurements 
at the upcoming LHC It is these more complicated processes where the power of our new 
formalism becomes increasingly pronounced. While the number and complexity of Dirac and 
color structures grows quickly for any exhaustive study of factorization with two or more 
final state jets, the application of our formalism is straightforward and in fact facilitates 
recycling the bulk of the work needed or already known in the literature for a particular 
observable, to be used for other observables of interest. 
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